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ABSTRACT 

We present an analysis of the 5-8 fim Spitzer-IRS spectra of a sample of 68 local 
Ultraluminous Infrared Galaxies (ULIRGs). Our diagnostic technique allows a clear 
separation of the active galactic nucleus (AGN) and starburst (SB) components in 
the observed mid-IR emission, and a simple analytic model provides a quantitative 
estimate of the AGN/starburst contribution to the bolometric luminosity. We show 
that AGNs are ~ 30 times brighter at 6 /im than starbursts with the same bolomet- 
ric luminosity so that even faint AGNs can be detected. Star formation events are 
confirmed as the dominant power source for extreme infrared activity since ~ 85% of 
ULIRG luminosity arises from the SB component. Nonetheless an AGN is present in 
the majority (46/68) of our sources. 
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1 INTRODUCTION 

Ultraluminous Infrared Galaxies (ULIRGs, Lm. > 1O 12 L0) 
are the local counterparts of the high-redshift objects dom- 
inating the cosmic background in the far-infrared and milli- 
metric bands. Unveiling the nature of their energy source is 
fundamental in order to understand the star formation his- 
tory and the obscured AGN activity in the distant Universe. 

Since their discovery, several infrared indicators have 
been proposed to determine whether the central engine in 
ULIRGs is an AGN or a starburst (SB). The presence of 
high- ionization lines in the mid-IR spectra of ULIRGs points 
to AGN activity, while intense PAH emission features are 
typical of starburst environments (Genzel et al. 1998; Lau- 
rent et al. 2000). Recently, the absorption feature of amor- 
phous silicate grains centered at 9.7 /an has also been used 
together with the PAH emission to assess the nature of the 
obscured power source (Spoon et al. 2007). An alternate way 
to disentangle AGNs and SBs in ULIRGs has been proposed 
by Risaliti et al. (2006, hereafter R06), based on the sepa- 
ration of the two continuum components in 3-4 /im spectra. 
This method has been successfully applied to a sample of 
~50 nearby ULIRGs (Risaliti, Imanishi & Sani 2007, sub- 
mitted) and provided an estimate of the average AGN/SB 
contribution to ULIRGs. The key reason for using the con- 
tinuum emission at A ~3-4 /xm as a diagnostic is the differ- 



ence between the 3-/tm to bolometric ratios in AGNs and 
starbursts (~two orders of magnitude larger in the former). 
This makes the detection of the AGN component possible 
even when the AGN is heavily obscured and/or bolometri- 
cally weak compared to the starburst. However the original 
prescription is limited by the low quality of the available L- 
band spectra of ULIRGs (R06, Imanishi et al. 2006), which 
makes the results on individual sources highly uncertain, 
except for the ~10-15 brightest objects. 

At present we have extended the analysis to the 5-8 /im 
spectral band, using the observations of the IRS instrument 
(Houck et al. 2004) onboard Spitzer. We disentangled the 
AGN and SB contributions to the observed 5-8 /im emission 
of ULIRGs by combining average spectral templates repre- 
senting the different properties of the two physical processes 
at work. The high quality of Spitzer-IRS data, in addition to 
the relatively low dispersion of the intrinsic continuum prop- 
erties of both AGNs and starbursts in this spectral range, 
allows a much more accurate determination of the AGN/SB 
components than possible at other wavelengths (e.g. X-rays) 
or with other diagnostic methods based on emission lines. 
In this paper we present our decomposition method, and 
discuss a simple analytical model providing a quantitative 
estimate of the AGN/SB contribution to the bolometric lu- 
minosity of each source. 



2 E. Nardini et al. 



2 OBSERVATIONS AND DATA REDUCTION 

In order to perform a detailed study of a representative sam- 
ple of ULIRGs in the local Universe, we selected 68 sources 
with z < 0.15 and a 60- /im flux density feo > 1 Jy- Most 
of the objects are taken from the IRAS ULIRG 1 Jy sam- 
ple (Kim & Sanders 1998), but a few more sources in the 
southern hemisphere have also been included. The flux limit 
at 60 /im ensures an unbiased selection with respect to the 
relative AGN/SB contributions. 

IRS observations were obtained within three differ- 
ent programs: PID 105 (PI J.R.Houck), PID 2306 (PI 
M.Imanishi), PID 3187 (PI S.Veilleux). The coadded im- 
ages provided by the Spitzer Science Center (after the treat- 
ment with pipeline version S13.0) have been background- 
subtracted by differencing the two observations in the nod- 
ding cycle. The spectra have been extracted and calibrated 
following the standard procedure for point-like sources with 
the package SPICE. The flux uncertainties have been esti- 
mated from source and background counts (in e~/s). Fi- 
nally, we performed a smooth connection between the Short- 
Low spectral orders, with no necessity of relative scaling. 

Out of the 68 spectra, we already published 48 in Iman- 
ishi et al. 2007. Six more spectra are shown in Armus et 
al. 2006, 2007. The remaining 14 spectra are analyzed here 
for the first time and will be fully presented in a forthcoming 
paper (Nardini et al. 2008, in prep.). 

3 THE 5-8 /im AGN/SB SEPARATION 

Despite the diversity of the global IRS spectra of pure AGNs 
and pure SBs, and the complexity of the physics involved, lit- 
tle dispersion is seen at wavelengths shortward of the 9.7 /im 
silicate feature. This makes possible the use of universal 
AGN/SB templates to reproduce the spectral properties of 
ULIRGs in the 5-8 fim interval. In the following we describe 
the templates adopted in our model. 

Starburst. The mid-IR spectral features of local SB 
galaxies show very little variations from one object to an- 
other in the 5-8 /xm wavelength range, while larger differ- 
ences are present in the ~9-30 fim band (Brandl et al. 2006, 
hereafter B06). In order to check if this is the case in the 
ULIRG luminosity range as well, we analyzed the sources 
in our sample estabilished to be starburst-dominated by 
multi-band studies. We did not find significant variations 
among the spectra, concluding that a fixed template can 
be used to represent the 5-8 /jm SB component in local 
ULIRGs. We built such template using the five brighest ob- 
jects among the pure SBs in our sample (IRAS 10190+1322, 
IRAS 12112+0305, IRAS 17208-0014, IRAS 20414-1651 
and IRAS 22491 — 1808), whose underlying continuum has 
been reproduced by a power law and normalized at 6 /xm 
before adding. Our SB template is shown in Fig[T] together 
with its dispersion in the whole IRS spectral band and the 
B06 template. The little spectral dispersion below 8 /an 
among SBs of different luminosity (to be compared with 
their large differences at longer wavelengths) is in itself an 
interesting result, which should be fully investigated through 
detailed emission and radiative transfer models. Concerning 
this we only notice that such a remarkable similarity can 
result from the spatial integration over a large number of 
individual star-forming regions. 
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Figure 1. Comparison between the SB template of B06 (dashed 
red line) and our template (solid blue line), constructed from the 
emission of the 5 brightest SB-dominated ULIRGs. The shaded 
area shows the 1-a rms dispersion in the 5 ULIRG spectra. The 
vertical long-dashed green lines enclose the fitting region. 



AGN. Our recent L-band analysis of bright ULIRGs 
shows that the intrinsic AGN emission is due to hot dust 
grains and the flux density is well described by a featureless 
power law with a fixed spectral index: /„ oc A . Here we 
adopt the same spectral shape up to 8 fim, in agreement with 
new Spitzer observations of a large sample of local type 1 
quasars (Netzer et al. 2007). 

An active nucleus is much more compact than a cir- 
cumnuclear starburst region. As a consequence, the near-IR 
radiation due to thin dust reprocessing can be itself strongly 
reddened because of a compact absorber along the line of 
sight. We therefore introduce an exponential attenuation 
factor e~ T< - x \ where the optical depth follows the conven- 
tional law t(A) oc A -1,75 (Draine 1989). A similar correc- 
tion is not needed in the SB template. We stress that this 
does not imply that the starburst spectrum is not affected 
by inner obscuration; the possible effects of this obscuration 
(which are clear at longer wavelengths, e.g. in the silicate ab- 
sorption features at 9.7 and ~18 /xm) are however already 
accounted for in the adopted observational template. 

Summarizing, the different contributions to the ob- 
served energy output of a ULIRG can be parametrized as 
follows: 

f: bs W = ft [(1 - a 6 K 6 + a 6 uT e- T(A) ] (1) 

where as is the AGN contribution to the 5-8 fim intrinsic 
flux density ft*, while uf and u a J> n are the SB and AGN 
templates normalized at 6 /im. Apart from the flux normal- 
ization, our model contains only two free parameters, i.e. ae 
and the optical depth to the AGN r(6 imi). They are both 
shown in Tab[T] 

Additional high-ionization emission lines and molecu- 
lar absorption features (due to ices and aliphatic hydrocar- 
bons), whenever present, were fitted by means of gaussian 
profiles except for the water ice absorption at ~6 fim, repro- 
duced with the laboratory profile from the Leiden database 
corresponding to pure H2O ice at 30 K. 
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Figure 2. Three representative examples of the typical 5—8 /an 
spectral shapes of ULIRGs. The differences among the spectra 
are entirely due to the different AGN contribution and its ob- 
scuration. Whenever the AGN is the dominant power source, as 
in Mrk 231, a strong continuum almost obliterates the PAH fea- 
tures. On the contrary, in Mrk 273 the AGN is fainter and the 
spectral outline of a starburst is clearly identified. A similar spec- 
trum is exhibited by IRAS 20551—4250, but the features are less 
prominent and the continuum is steeper: this source harbours an 
obscured AGN. In each panel, in addition to the data (green filled 
circles) and their best fits (black thin line), we have included the 
reddened AGN (red dot-dashed line) and starburst (blue dotted 
line) components. 



FigE] shows the spectral decomposition of three rep- 
resentative ULIRGs. In spite of the great diversity of the 
observed spectra, our simple model provides a good fit of 
each spectrum in the sample: the residuals from the best 
fits are in all sources smaller than 10% at all wavelengths 
(though the fits are not formally acceptable in a statistical 
sense, with a reduced x 2 ^2, due to the small error bars 
and the remaining unfitted minor features). In particular, 
both the PAH emission and the continuum are always well 
reproduced: this implies that the large variations in the 5- 
8 fj,m spectral shape of ULIRGs are entirely due to the AGN 
contribution and its obscuration. A detailed analysis of the 
results for each source and a physical interpretation of pe- 
culiar cases are the subject of a forthcoming paper (Nardini 
et al. 2008, in prep.). 



4 AGN/SB BOLOMETRIC CONTRIBUTIONS 

The large difference between the 5-8 fj,m to bolometric ratios 
in AGNs and SBs implies that this ratio is itself an indicator 



of AGN activity, and can be used (a) to test the consistency 
of our decomposition method and (b) to estimate the relative 
AGN and SB contributions to the bolometric luminosity of 
our sample. We define the 6-/xm to bolometric ratio as: 



R = 



pint 

"ah 
Ftr 



(2) 



where Fin is the total infrared flux, estimated as in Sanders 
& Mirabel (1996). Since the integrated infrared luminosity of 
ULIRGs coincides almost exactly with their bolometric lu- 
minosity, R is a fair approximation to the fraction of the to- 
tal energy output that is intrinsically emitted in the 5-8 /im 
range. Reminding that the intrinsic AGN/SB contributions 
are aefe nt and (1 — ae)fe nt respectively, and decomposing 
Fir as Fj^ n + Fjn, a simple connection between R and ae 
is brought out: 



R 



R a 9 n R sb 



aaR sb + (1 - a 6 )R a 



(3) 



provided that R agn and R sb , the equivalents of R for pure 
(unobscured) AGNs and pure SBs, are defined as in Eq0 
We fitted the theoretical R(aa) relation (Eq(3J to our data 
considering R agn and R sb as free parameters, and found: 



logR agn = -0.49±C.i3 and 



logi? s6 = 



qo+o.03 



(4) 



We note that R a s n turns out to be somewhat higher than tra- 
ditional estimates based on AGN spectral energy distribu- 
tions: for example, we derive logi?" 9 " ~ —0.6 from the SED 
of Elvis ct al. (1994). This suggests that the local quasars 
(mostly PG quasars) used to build the mentioned SED can 
be contaminated to some extent by a starburst contribution, 
in agreement with recent studies (Netzer et al. 2007). 

According to Eq(?l AGNs are ~30 times more luminous 
at 6 (im than starbursts with the same bolometric lumi- 
nosity. We are now able to quantify the AGN contribution 
(ctboi = Fj^ /Fm) to the total infrared luminosity of each 
source: 



Cibol = 



a 6 + {R agn /R sb )(l~a 6 )- 



(5) 



where R a<3n / R sb ~ 28. The values of a bol are listed in 
Tab[T] Our estimates for the ~15 brightest sources are in 
good agreement with those of R06 and with the Genzel et 
al. (1998) and Laurent et al. (2000) mid-IR diagnostic di- 
agrams. Considering the whole sample, our results can be 
compared with the optical classification and with L-band 
and hard X-ray studies, when available. A substantial agree- 
ment is obtained in all cases. It is worth noticing that the 
optical classification alone gives incomplete information: all 
the sources classified as Seyferts show clear traces of AGN 
activity, but 7 out of 8 among the ULIRGs with abol > 0.25 
and r > 1 are indeed classified as LINERs or H II regions. 
LINERs are again confirmed to be rather heterogeneous with 
respect to the nature of their energy source. Such ambigui- 
ties can be solved by applying our diagnostic. 

By inverting Eq[5]the relation between R and «(,„; takes 
the neat form R = a ho iR agn + (1 - a bo i)R sh ', and is plotted 
in Fig(3^. As a final check we have computed for each source 
the following quantities: 



R" 



\ aboiFiR 



and R 



fe(l 



(1 - a bo i)F m 



(6) 
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Table 1. Spectral parameters for the 68 sources in our sample, qq: AGN contribution to the intrinsic continuum emission at 6 fira 
(in percent), r: Optical depth of the AGN component at 6 /an (we assume r = for the sources with no detected AGN). ctbol'- AGN 
contribution to the bolometric luminosity (in percent). The errors in a 00 i are due to the statistical uncertainty both in the flux amplitude 
of the AGN/SB components and in the ratios R agn and R sb . The systematic effects are discussed in the text. O/X/L: SB/ AGN /LINER 
classification based on optical, X-ray and L-band spectroscopy. A: AGN, L: LINER, A*: AGN, tentative detection. References: : Vcillcux 
et al. 1999, 2 : Veilleux et al. 1995, 3 : Due et al. 1997, 4 : Iwasawa et al. 2005, 5 : Severgnini et al. 2001, 6 : Franceschini et al. 2003, 7 : Balestra 
et al. 2005, 8 : Vignati et al. 1999, 9 : Imanishi et al. 2003, 10 : Imanishi et al. 2006, 11 : Risaliti et al. 2006, 12 : Sani et al. 2007. + : Sources 
with simultaneous ajoj > 0.2 and r > 1. 
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The results are shown in Fig(3]3 and prove that our de- 
composition method is reliable in estimating the AGN/SB 
contributions both to the 6-/xm and to the bolometric lu- 
minosity of local ULIRGs. In fact, the estimated 6-/xm to 
bolometric ratios of the SB component in the composite 
sources (which are located at the bottom right of the plot 
and are heavily dependent on our modeling) are fully consis- 
tent, within the errors, with the ratios of the pure starbursts 
(located at the bottom left and directly computed from the 
measured 6-/im and IRAS fluxes). This success is promis- 
ing in anticipation of a forthcoming study about the role 
of black hole accretion and star formation in the intense 
infrared activity characterising the distant galaxies. 

However, it is important to keep in mind the main lim- 
itations of our approach: 

1) The narrow wavelength range used in this work simplifies 
the decomposition analysis, but prevents us from a complete 
study of the dust composition, density and geometrical dis- 



tribution. These elements strongly affect the overall mid-IR 
emission longward of the silicate absorption feature, and can 
be investigated only through an analysis of the whole IRS 
spectrum. 

2) While the 5-8 /im templates seem to have little dispersion 
(as discussed in detail in Section |3}, the spread in the 6-[im 
to bolometric ratios R a3n and R sb can be much higher, mak- 
ing our estimates of the AGN/SB bolometric fractions more 
uncertain than those in the 5-8 /im band. The uncertainties 
in ctboi reported in Tab[T] are obtained assuming the mean 
ratios, with the errors on the mean given in Eq|3] However, 
the dispersion around the best fit in Fig(3^ is significantly 
larger. We therefore consider this dispersion (0.3 dex, con- 
stant at all values of Q„ i) as the actual uncertainty in the 
bolometric ratios of the individual sources. The numerical 
results on the individuale sources are anyway precise enough 
to estabilish which is the dominant source of the observed 
luminosity. 
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Figure 3. (a) Ratio R between absorption-corrected 6-/mi lu- 
minosity and bolometric luminosity, versus the AGN bolometric 
contribution ajoi- The error bars of R are due to the uncertainties 
in the total infrared flux Fjh and in the intrinsic AGN fraction, 
ae- The solid line is the best fit of the R— a% relation from Eq|3] 
(plotted as a function of using Eq[5}. (b) Same as above, 
with the AGN and SB components plotted separately. 



Overall, if we consider as a confidence limit the value 
aboi = 0.01 (i.e. the dispersion around our starburst tem- 
plate), an AGN is present in 46 of the 68 ULIRGs in our 
sample (including several of those optically classified as H II 
regions), but it is significant (atoi > 0.25) only in ~30% of 
the cases. The SB process is responsible for almost 90% of 
the observed infrared luminosity of ULIRGs, with no signif- 
icant (i.e. > 5%) bias due to the sample selection. A simi- 
lar fraction holds for the subsample of 34 sources optically 
classified as LINERs. Our analysis is also consistent with 
the findings about the nature of high-redshift infrared-bright 
galaxies detected in 24 /im Spitzer MIPS surveys. IRS spec- 
troscopy shows that they are mostly z ~1— 3 galaxies, with 
an apparent bias toward AGN-dominated sources (Houck et 
al. 2005). This is in agreement with the ~30 times higher 
AGN relative emission in the 5-8 fj,m rest-frame wavelength 
range we have pointed out. 



spectra of 68 local ULIRGs, observed with the Spitzer Space 
Telescope. We have been able to detect an AGN in more 
than 60% of our sources, and estimate its contribution to 
the bolometric luminosity. In a statistical sense, we con- 
firm that local Ultraluminous Infrared Galaxies are powered 
for ~85% by intense star formation and for the remaining 
~15% by AGN activity. Our method proves to be successful 
in unveiling an intrinsically faint or obscured AGN inside a 
ULIRG. In this context we also put on a sound basis our 
initial assumption that the wavelength interval 5-8 /im is 
an appropriate spectral range in order to search for AGNs: 
an AGN turns out to be approximately 30 times more lu- 
minous at 6 fim than a starburst with the same bolometric 
luminosity. 
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5 CONCLUSIONS 



The use of average templates for AGN and SB emission has 
allowed us to disentangle the two components in the 5-8 /im 



